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O ' Abstract 
(N 

5 ' The decay rates and spectroscopy of the D and Ds mesons are computed in a nonrela- 

tivistic phenomenological quark-antiquark potential of the type V{r) = — |^ + Ar'^ with 

different choices of v. Numerical method to solve the Schrodinger equation has been used 

to obtain the spectroscopy of qQ mesons. The spin hyperfine, spin-orbit and tensor com- 

i-^ ' ponents of the one gluon exchange interactions are employed to compute the spectroscopy 

Oh' 

I ' of the excited S states, low lying P-waves and D-waves. The numerically obtained radial 

D-i' 

QJ solutions are employed to obtain the decay constant and leptonic decay widths. It has been 

observed that predictions of the spectroscopy and the decay widths are consistent with 
other model predictions as well as with the known experimental values. 
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'. 1 Introduction 

oo 

o 
a^ 
o 



Spectroscopy of hadrons containing heavy flavours has attracted considerable interest in recent 
years due to many experimental facilities such as the BES at the Beijing Electron Positron 
CoUider (BEPC), E835 at Fermilab, and CLEO at the Cornell Electron Storage Ring (CESR) 
etc., worldover. They have been able to collect huge data samples in the heavy flavour sector. 
Where as B-meson factories, BaBar at PEP-II and Belle at KEKB are working on the observation 
of new and possibly exotic hadronic states. All these experiments are capable of discovering new 
hadrons, new production mechanisms, new decays and transitions and in general will be providing 
high precision data sample with better stastics and higher confidence level. After having played 
a major role in the foundation of QCD, heavy hadron spectroscopy has witnessed in the last few 
years a renewal of interest led by the many new data coming from the B factories, CLEO and 
the Tevatron and by the progress made in the theoretical methods. The remarkable progress at 
the experimental side, with various high energy machines such as BaBar, BELLE, B-factories, 
Tevatron, ARGUS collaborations, CLEO, CDF, D0 etc., for the study of hadrons has opened 
up new challenges in the theoretical understanding of light-heavy flavour hadrons. The existing 
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results on excited heavy-light mesons are therefore partially inconclusive, and even contradictory 
in several cases. The predictions of masses of heavy-light system for ground state as well as excited 
state are few from the theory [HEjlSllllElEllT]. In the open charm sector, the observation of 
a charm-strange state, the D*j(2317) state [8j by BaBar Collaboration. It was confirmed by 
CLEO Collaboration at the Cornell Electron Storage Ring p] and also by Belle Collaboration 
at KEK [To]. Besides, BaBar had also pointed out to the existence of another charm-strange 
meson, the Dsj(2460) [8]. This resonance was measured by CLEO [9J and confirmed by Belle 
[10] . Belle results [lO] are consistent with the spin-parity assignments of = 0^ for the 
£1*^(2317) and = 1+ for the L),j(2460). Thus, these states are well established and confirmed 
independently by different experiments. They present unexpected properties, quite different 
from those predicted by quark potential models. If they would correspond to standard P— wave 
mesons made of a charm quark and a strange antiquark their masses would be larger [TT], around 
2.48 GeV for the D*j{2317) and 2.55 GeV for the Dsj{24m). They would be therefore above the 
DK and D*K thresholds, respectively with being broad resonances. However the states observed 
by BaBar and CLEO are very narrow, T < 4.6 MeV for the D*j{2317) and T < 5.5 MeV for the 
Dsj{2m). 

In near future, even larger data samples are expected from the BES-III upgraded experiments, 
while the B factories and the Fermilab Tevatron will continue to supply valuable data for few 
more years. Later on, the LHC experiments at CERN, Panda at GSI etc., will be accumulating 
large data sets which will offer greater opportunities and challenges particularly in the field of 
heavy flavour physics [T2] . 

At the hadronic scale the nonperturbative effects connected with complicated structure of QCD 
vacuum necessarily play an important role. But our limited knowledge about the nonperturbative 
QCD leads to a theoretical uncertainty in the quark- antiquark potential at large and intermediate 
distances [I3]. So a successful theoretical model can provide important information about the 
quark- antiquark interactions and the behavior of QCD at the hadronic scale. Though there 
exist many potential models with relativistic and nonrelativistic considerations employed to 
study the hadron properties based on its quark structure [HI [151 [HI [13 [HI [HI [SHI [211 [221 
[23l [24l [25l [26l [271 [28] . the most commonly used potential is the coulomb plus linear power 
potential, V{r) = + err, with the string tension a [221 [20] • However, for the higher excited 

mesonic states it is argued that the string tension a must depend on the QQ separation [ST] [32]. 
This corresponds to flattening of the confinement potential at larger r (r > Ifm). More over 
the analysis based on Regge trajectories for meson states suggests the confinement part of the 
potential to have the power | instead of 1 [331 [31] • This has prompted us to choose a power 
form for the confining part of the interquark potential and study the properties of heavy flavour 
systems by varying the power of the confinement part of the interquark potential different from 
1.0. 
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Apart from the spectroscopic predictions of higher orbital states, other problems associated with 
the phenomenological models employed for the properties of mesons are the right predictions 
of their decay properties. For better predictions of the decay widths, many models need to 
incorporate additional terms such as the radiative contributions, higher order QCD corrections 
etc., to the conventional decay formula |T1[ |15l [35], [361 [37] . 

The decay widths can provide an account of the compactness of the meson system in terms of the 
radial wave function which is an useful information complementary to spectroscopy [3H]. Other 
unresolved issues are related to the hyperfine and fine structure splitting of the mesonic states and 
their intricate dependence with the constituent quark masses and the running strong coupling 
constant. Thus, in this paper we make an attempt to study properties like mass spectrum, decay 
constants and other decay properties of the open charm mesons {D,Ds). We investigate the 
heavy- light mass spectra of D{cq) and Ds{cs) mesons in the frame work of the nonrelativistic 
CPPjy potential model. In the present study, we consider different choices of the potential power 
index u in the range 0.1 < < 2.0. 



2 Nonrelativistic Treatment for Heavy Flavour Mesons 
using CPP^y 



In general, properties of heavy flavour mesons have been studied based on potential models 
in the frame work of relativistic as well as nonrelativistic quantum mechanics. In the limit of 
heavy quark mass tuq oo, heavy meson properties are governed by the dynamics of the light 
quark. As such, these states become hydrogen like atoms of hadron physics. Moreover, both the 
non-relativistic predictions are in fair agreements with each other as well as with the available 
experimental and lattice results. Hence, for the present study of charm meson bound states, we 
consider a nonrelativistic Hamiltonian given by [Til [T5l [391 [lOl [Hj [12] 

^ = + 2^ + '''<'■'• <1> 

here, 

M = mi + 1712, and Mi = , (2) 

nil + 

The relative momentum of each quark is represented by p and V{r) is the quark- antiquark 
potential. Nonrelativistically, this interaction potential consists of a central term V^(r) and a 
spin dependent part Vsoi.^)- The central part Vc{r) is expressed in terms of a vector (Coulomb) 
plus a scalar (confining) part given by 

4 (T 

V,{r) = Vy + Vs = --- + AT^ (3) 

3 r 
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as the static quark-antiquark interaction potential [H]. 




2.1 Spin-dependent Forces 

In general, the quark-antiquark bound states are represented by n'^^'^^Lj, identified with the 
J^^ values, with J = L + S, S = Sq + Sq, parity P = (—1)-^+^ and the charge conjugation 
C = (— l)'^"^'^ with {n,L) being the radial quantum numbers. So the S'-wave (L = 0) bound 
states are represented by J'^'" = 0~+ and 1 respectively. The P-wave (L = 1) states are 
represented by J^'" = with L = 1 and S = while J^'" = 0"'"+, l^"*" and 2+"*" correspond 
to L = 1 and S = 1 respectively. Thus, the spin-spin interaction among the constituent quarks 
provides the mass splitting of J = ^ and 1 states, while the spin-orbit interaction provides 
the mass splitting of J^*" = 1+^ and 2"''+ states. The J^'-' = state with L = 1 and 

S = represents the spin average mass of the P-state as its spin-orbit contribution becomes 
zero, while the two J = singlet and the J = l^"*" of the triplet P-states can form a mixed 
state. The D-wave (L = 2) states are represented by J^'-' = 2 ^ with L = 2 and S = while 
J^'^ = 3 , 2 and 1 correspond to L = 2 and S = 1 respectively. The P-wave (L = 3) 
states are represented by J^'" = 3+^ with L = 3 and S = while J^'~' = 4++, 3+^ and 2+^^ 
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correspond to L = 3 and S = 1 respectively. 



For computing the mass difference between different degenerate meson states, we consider the 
spin dependent part of the usual one gluon exchange potential (OGEP) given by [271 SSI HH 113 
16] . Accordingly, the spin-dependent part, Vsni^) contains three components of the interaction 
terms, such as the spin-spin, the spin-orbit and the tensor part given by 



VsD{r) = Vss{r] 
Vrir) 



5(5+1) 



+ VLsir) {L-S) + 
3{S-r){S-r) 



(4) 



The spin-orbit term containing VLS'(r) and the tensor term containing Vrir) describe the fine 
structure of the meson states, while the spin-spin term containing V5s(r) proportional to 2(s^ • 
s'q) = S{S + 1) — I gives the spin singlet-triplet hyperfine splitting. 

The coefficient of these spin-dependent terms of EqnJH can be written in terms of the vector 
and scalar parts of the static potential, V^(r) as 



T/ r ^ 1 f.dVv dVs\ 

VT{r) = 3-— — (6) 

mim2 \ dr'^ r dr J 

Vssir) = :-^VVv = i^5(^)(rl (7) 

3 7711777,2 9 7771777,2 

The present study with the choices of u in the range 0.1 < < 2.0, is an attempt to know the 
predictability of the hadron spectroscopy with a chosen value of mass parameters (7771,7772) and 
confinements strength represented by the potential parameter A. The running strong coupling 
constant appeared in the potential V{r) in turn is related to the quark mass parameter as 

4:71 

> = (11 - In,) 1„(,VA^) 

Where, nj is the number of flavors, fi is renormalization scale related to the constituent quark 
masses as /7 = 2mim2/{mi +1712) and A is the QCD scale which is taken as 0.150 GeV by fixing 
as = 0.118 at the Z-boson mass (91 GeV)^. 

The potential parameter, A of Eqnj3] is similar to the string strength a of the Cornell potential. 
The different choices of u here then correspond to different potential forms. So, the potential 
parameter A expressed in GeV^^ can be different for each choices of u. The model potential 
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Table 1: Square of the radial wave functions at the origin (|i?„(0)p (in GeV'^)) of Qq systems in 
CPP,. 



D D, 



Model 



1 






IS 




2S 




3S 




IS 




2S 




3S 


0. 


,1 


0. 


,050 


0, 


.011 


0, 


.005 


0, 


.063 





.014 





.007 


0. 


,3 


0. 


,083 


0, 


.028 


0, 


.016 


0, 


.104 





.034 





.020 


0. 


,5 


0. 


,111 


0, 


.047 


0, 


.031 


0, 


.139 





.058 





.038 


0. 


,7 


0. 


,135 


0, 


.068 


0, 


.049 


0, 


.168 





.084 





.061 


0, 


,8 


0, 


,145 


0, 


.079 


0, 


.060 


0, 


.181 





.098 





.074 


0. 


,9 


0, 


,155 


0, 


.091 


0, 


.071 


0, 


.193 





.112 





.088 


1. 


,0 


0. 


,164 


0, 


.102 


0, 


.083 


0, 


.204 





.127 





.103 


1. 


,1 


0. 


,172 


0, 


.114 


0, 


.096 


0, 


.214 





.142 





.118 


1. 


,3 


0. 


,186 


0, 


.139 


0, 


.123 


0, 


.232 





.172 





.153 


1. 


,5 


0. 


,199 


0, 


.164 


0, 


.154 


0, 


.247 





.203 





.190 


2. 


,0 


0. 


,222 


0, 


.227 


0, 


.237 


0, 


.276 





.280 





.292 



parameter A and the mass parameter of the quark/antiquark (mi, 1712) are fixed using the known 
ground state center of weight (spin average) mass and the hyperfine splitting (Ms 5^ — Mig^) of 
D and Dg systems respectively. The spin average mass for the ground state is computed for the 
different choices of v in the range, 0.1 < z/ < 2.0. The spin average or the center of weight mass. 
Mew is calculated from the known experimental/theoretical values of the pseudoscalar (J = 0) 
and vector ( J = 1) meson states as 

E(2J + 1) M^j 
"'-^'^ = 'e(2J+i) 

J 



2.2 Spectra of Heavy - Light Flavour {Qq) Mesons 

The spectra of the heavy-light mesons are calculated using nonrelativistic hamiltonian as given 
by Eqn{T], where mi = uiq and 1712 = rrig. The spin average masses of D* — D and the D* — Dg 
mesons are computed using the experimental values of Md = 1.869 GeV, M^. =2.010 GeV, 
Md, =1.968 GeV, M^j = 2.112 GeV respectively [17|. 

We employ the numerical approach as given by [H] to find the eigen values and radial wave 
functions of the respective Schrodinger equation. The potential parameter A, is made to vary 
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Table 2: The derivative of orbitally excited radial wave functions at the origin (|i?^(0)| in 

GeV (i+^)) of Qq systems in CPP^. 

D D, 



Model 



[ 






IP 




2P 




ID 




IF 




IP 




2P 




ID 




IF 


0, 


.1 


0, 


.0110 


0, 


,0083 


0, 


.0008 





.0001 


0, 


.0130 


0, 


.0097 


0, 


.0010 


0, 


,0001 


0, 


.3 


0, 


.0253 


0, 


,0210 


0, 


.0034 





.0006 


0, 


.0299 


0, 


.0253 


0, 


.0044 


0, 


,0008 





.5 


0, 


.0392 


0, 


,0358 


0, 


.0076 





.0018 


0, 


.0474 


0, 


.0428 


0, 


.0097 


0, 


,0024 


0, 


.7 


0, 


.0551 


0, 


,0520 


0, 


.0129 





.0040 


0, 


.0634 


0, 


.0619 


0, 


.0165 


0, 


,0055 


0, 


.8 


0, 


.0603 


0, 


,0604 


0, 


.0159 





.0055 


0, 


.0711 


0, 


.0722 


0, 


.0203 


0, 


,0076 


0, 


.9 


0, 


.0666 


0, 


,0692 


0, 


.0192 





.0073 


0, 


.0793 


0, 


.0826 


0, 


.0245 


0, 


,0103 


1, 


.0 


0, 


.0732 


0, 


,0780 


0, 


.0226 





.0093 


0, 


.0881 


0, 


.0931 


0, 


.0288 


0, 


,0126 


1, 


.1 


0, 


.0790 


0, 


,0870 


0, 


.0261 





.0115 


0, 


.0934 


0, 


.1029 


0, 


.0335 


0, 


,0157 


1, 


.3 


0, 


.0916 


0, 


,1043 


0, 


.0334 





.0166 


0, 


.1088 


0, 


.1243 


0, 


.0426 


0, 


,0227 


1, 


.5 


0, 


.1022 


0, 


,1221 


0, 


.0408 





.0225 


0, 


.1189 


0, 


.1459 


0, 


.0521 


0, 


,0307 


2, 


.0 


0, 


.1234 


0, 


,1653 


0, 


.0593 





.0397 


0, 


.1457 


0, 


.1969 


0, 


.0756 


0, 


,0542 



with z/, keeping the quark mass parameter fixed for each choices of Qq system. It is observed that 
the hyperfine splitting of the 1^5*1 and 1^5*0 states are very sensitive to the choices of uiq and A. 
The most suitable values of the quark mass parameter are found to be = 1.28 GeV, nifi = 0.35 
GeV rris = 0.500 GeV for our present study. The corresponding A values obtained from the 15* 
fitting and are plotted in Fig{T] against the potential index u of the D and systems. Just 
like the string tension a{r) of the Cornell potential was made to vary for excited states [311 132], 
we allow A to vary mildly with radial quantum number n as A ^ for computing the spin 

(n)3 

independent masses of the orbital excited (nL) states. The variation in A can be justified by 
similar arguments for the changes in with the average kinetic energy. Here, as the system 
get excited, the average kinetic energy increases and hence the potential strength (the spring 
tension) reduces. With this mild state dependence on the potential parameter A, we obtain the 
spin average masses of the orbital excited states closer to the experimentally known D and 
systems. The computed values of the radial wave function at the origin |i?^£(0)| for different 
states are listed in Tabled] (nS"— states) and Table [2] (IP, 2P, ID, IF— states) for all the cq {Q e c 
and q e u/d,s) combinations. Using the spin dependent potential given by EqnJH we compute 
the masses of the different low lying states of cq and are listed in Table [3] and H] in the 

case of D and Dg mesons respectively. The available experimental values as well as other model 
predictions are also listed for comparison. 
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Table 3: Mass spectra (in GeV) of D meson. 



State 










Potential index u 












Expt. 

[47] 


RQM 

m 


RQM 

m 


BSU 

[49] 


0.1 


0.7 


0, 


.8 


0.9 


1.0 


1.1 


1.3 


1.5 


1.7 


2.0 




1.985 


2.007 


2, 


.010 


2.013 


2.015 


2.018 


2. 


.021 


2.025 


2.028 


2.031 


2.010 


2.009 


2.005 


2. 


,006 




1.932 


1.864 


1, 


.855 


1.848 


1.841 


1.834 


1, 


.823 


1.813 


1.805 


1.794 


1.869 


1.875 


1.868 


1. 


,874 




2.070 


2.268 


2, 


.294 


2.319 


2.342 


2.364 


2, 


.404 


2.440 


2.472 


2.514 


2.460 


2.459 


2.460 


2. 


,477 




2.072 


2.282 


2, 


.310 


2.336 


2.361 


2.384 


2, 


.426 


2.465 


2.498 


2.542 




2.414 


2.417 


2. 


,407 




2.068 


2.261 


2, 


.287 


2.312 


2.335 


2.357 


2, 


.398 


2.434 


2.467 


2.510 




2.438 


2.490 


2. 


,341 


J- -< 1 




9 91 fi 


9 


.zoo 


2.253 


2.269 


2.285 


2, 


.312 


9 '?'?7 
z.oo ( 


z.ooo 


9 

z.ooo 




9 p;ni 

Z.OUi 


9 '^77 


2. 


,389 




2.011 


2.303 


2, 


.350 


2.398 


2.445 


2.491 


2, 


.582 


2.668 


2.751 


2.868 




2.629 


2.692 


2, 


,601 


2'So 


1.998 


2.226 


2, 


.261 


2.296 


2.329 


2.362 


2, 


.425 


2.483 


2.538 


2.612 




2.579 


2.589 


2. 


,540 


I'D, 


2.106 


2.456 


2, 


.508 


2.558 


2.605 


2.651 


2, 


.736 


2.816 


2.887 


2.984 






2.799 


2. 


,688 




2.106 


2.454 


2, 


.504 


2.552 


2.597 


2.639 


2, 


.717 


2.788 


2.851 


2.933 






2.775 


2. 


,727 




2.105 


2.459 


2, 


.512 


2.564 


2.613 


2.661 


2, 


.751 


2.836 


2.913 


3.019 






2.833 


2. 


,750 




2.104 


2.455 


2, 


.509 


2.561 


2.612 


2.662 


2, 


.756 


2.846 


2.930 


3.046 






2.795 


2. 


,689 




2.037 


2.432 


2, 


.499 


2.565 


2.631 


2.696 


2, 


.824 


2.949 


3.069 


3.241 






3.035 


2. 


,860 




2.038 


2.443 


2, 


.511 


2.580 


2.648 


2.715 


2, 


.847 


2.976 


3.100 


3.279 






2.995 


2. 


,802 


23Po 


2.036 


2.427 


2, 


.493 


2.559 


2.624 


2.689 


2, 


.816 


2.940 


3.059 


3.231 






3.045 


2. 


,758 


ol p 


z.uoo 




9 




2.512 


2.570 


2.627 


2, 


.736 




9 QAA 


O.Uo ( 






9 QzLQ 


2. 


,792 




2.005 


2.468 


2, 


.553 


2.640 


2.727 


2.816 


2, 


.994 


3.172 


3.348 


3.607 






3.226 


2, 


,947 


3 5*0 


1.999 


2.413 


2, 


.486 


2.560 


2.634 


2.708 


2, 


.855 


2.999 


3.139 


3.340 






3.141 


2. 


,904 




2.127 


2.607 


2, 


.683 


2.757 


2.829 


2.900 


3, 


.033 


3.159 


3.276 


3.436 






3.091 








2.127 


2.599 


2, 


.671 


2.741 


2.807 


2.871 


2, 


.989 


3.096 


3.192 


3.318 






3.074 








2.127 


2.611 


2, 


.688 


2.764 


2.839 


2.911 


3, 


.051 


3.184 


3.308 


3.481 






3.123 








2.127 


2.617 


2, 


.698 


2.778 


2.857 


2.936 


3, 


.089 


3.238 


3.382 


3.587 






3.101 







Relativistic Quark Model (RQM), Blankenbecler- Suger Equation (BSU). 



Table 4: Mass spectra (in GeV) of Ds meson. 



State Potential index v Expt. RQM RQM BSU 





0.1 


0.7 


0, 


.8 


0.9 


1.0 


1.1 


1.3 


1.5 


1.7 


2.0 


133 


[26] 


12J 




E9] 




2, 


.086 


2.102 


2, 


.104 


2.106 


2, 


.108 


2.109 


2.112 


2.114 


2.116 


2.119 


2.112 


2.111 


2.113 


2. 


,108 




2, 


.047 


1.998 


1, 


.992 


1.987 


1, 


.982 


1.977 


1.969 


1.962 


1.956 


1.948 


1.968 


1.981 


1.965 


1. 


,975 




2, 


.165 


2.348 


2, 


.372 


2.394 


2. 


.416 


2.436 


2.474 


2.506 


2.535 


2.573 


2.572 


2.560 


2.581 


2, 


,586 




2, 


.162 


2.332 


2, 


.355 


2.376 


2. 


.397 


2.416 


2.452 


2.484 


2.512 


2.549 


2.535 


2.515 


2.535 


2. 


,522 




2, 


.157 


2.300 


2, 


.317 


2.334 


2, 


.350 


2.364 


2.391 


2.414 


2.433 


2.459 


2.317 


2.569 


2.487 


2. 


,455 


1 Pi 


2, 


.163 


2.337 


2, 


.360 


2.382 


2, 


.402 


2.422 


2.457 


2.488 


2.516 


2.552 


2.460 


2.508 


2.605 


2, 


,502 




2, 


.110 


2.355 


2, 


.395 


2.434 


2, 


.474 


2.513 


2.588 


2.661 


2.730 


2.828 




2.716 


2.806 


2. 


,722 




2, 


.101 


2.303 


2, 


.334 


2.365 


2, 


.395 


2.425 


2.482 


2.535 


2.585 


2.654 




2.670 


2.700 


2, 


,659 




2, 


.195 


2.502 


2, 


.545 


2.587 


2, 


.627 


2.666 


2.736 


2.799 


2.880 


2.965 






2.925 


2. 


,857 




2, 


.195 


2.505 


2, 


.551 


2.595 


2, 


.639 


2.681 


2.760 


2.832 


2.855 


2.929 






2.900 


2. 


,856 




2, 


.194 


2.502 


2, 


.548 


2.594 


2, 


.638 


2.681 


2.763 


2.839 


2.899 


2.989 






2.913 


2. 


,845 




2, 


.195 


2.503 


2, 


.548 


2.591 


2, 


.633 


2.674 


2.749 


2.818 


2.909 


3.007 






2.953 


2. 


,838 


23P2 


2, 


.136 


2.489 


2, 


.549 


2.608 


2, 


.668 


2.727 


2.843 


2.954 


3.040 


3.190 






3.157 


2. 


,988 




2, 


.134 


2.478 


2, 


.535 


2.593 


2, 


.651 


2.708 


2.820 


2.928 


3.062 


3.217 






3.114 


2, 


,942 


ri3 D 

2^Po 


2, 


.132 


2.455 


2, 


.507 


2.560 


2, 


.612 


2.663 


2.763 


2.858 


3.032 


3.182 






3.067 


2. 


,901 


2^Pi 


2, 


.135 


2.482 


2, 


.540 


2.598 


2, 


.656 


2.713 


2.826 


2.935 


2.949 


3.078 






3.165 


2. 


,928 




2, 


.105 


2.497 


2, 


.568 


2.641 


2, 


.715 


2.790 


2.940 


3.090 


3.238 


3.456 






3.345 


3. 


,087 


3 5*0 


2, 


.101 


2.459 


2, 


.523 


2.587 


2, 


.652 


2.717 


2.846 


2.973 


3.096 


3.275 






3.259 


3. 


,044 


13P4 


2, 


.214 


2.630 


2, 


.693 


2.755 


2, 


.814 


2.872 


2.979 


3.075 


3.221 


3.361 






3.220 








2, 


.214 


2.638 


2, 


.705 


2.771 


2, 


.836 


2.900 


3.022 


3.136 


3.163 


3.278 






3.224 






1=^P2 


2, 


.214 


2.642 


2, 


.711 


2.780 


2, 


.849 


2.917 


3.048 


3.175 


3.244 


3.392 






3.247 






ilPa 


2, 


.214 


2.635 


2, 


.702 


2.766 


2, 


.830 


2.892 


3.010 


3.119 


3.295 


3.466 






3.203 







Relativistic Quark Model (RQM), Blankenbecler- Suger Equation (BSU). 



3 The Decay constants of the charm flavored mesons 



The decay constants of mesons are important parameters in the study of leptonic or non-leptonic 
weak decay processes. The decay constants of pseudoscalar (/p) and vector (fy) states are 
obtained by parameterizing the matrix elements of weak current between the corresponding 
mesons and the vacuum as [50| 

{0\qri5c\P,{k))=zfpk^ (10) 

{0\qrc\V{k,e)) = fvMve'^ (11) 

where k is the meson momentum, e'^ and My are the polarization vector and mass of the vector 
meson. 

In the relativistic quark model, the decay constant can be expressed through the meson wave 
function $p/y(p) in the momentum space [26j. 



fp/v 



1/2 / p / N , _ \ 1/2 



12 f (Pp ( Ec{p) +mc\ ( Eg{p) + m, 



Mp/v J (27r)3 V 2E,{p) / V 2Eq{p) 



p2 



p2 

with Xp = —1 and Ay = —1/3. In the nonrelativistic limit — - << 1.0, this expression reduces to 
the well known relation between fp/v and the ground state wave function at the origin Rp/v{0) 
the Van-Royen-Weisskopf formula [^I]. Though most of the models predict the meson mass 
spectrum successfully, there exist wide range of predictions of their decay constants. For example, 
the ratio was predicted to be > 1 in most of the nonrelativistic cases, as mp < my and their 
wave function at the origin has assumed to be as -Rp(O) ~ -Ry(O) [52j. The ratio computed in the 
relativistic models [53] have predicted < 1, particularly in the QQ sector, but > 1 in the 
heavy-light flavour sector. The disparity of the predictions of these decay constants play decisive 
role in the precision measurements of the weak decay parameters. The value of the radial wave 
function (Rp) for 0^+ and (Ry) for 1 states would be different due to their spin dependent 
hyperfine interaction. The spin hyperfine interaction of the heavy flavour mesons are small and 
this can cause a small shift in the value of the wave function at the origin. Though, many 
models neglect this difference between (Rp) and (Rv), we consider this correction by making an 
ansatz that the -Rp/y(0) are related to the value of the radial wave function at the origin, Rn{0) 
according to the same way their masses are related. Thus, by considering 



MnP/V — Mn^c'W 



1 + {SF)p/y-^^''^' 



n,CW 



(13) 
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Table 5: The decay constants /p/y in MeV of D and Dg systems (The bracketed quantities are 
with QCD correction). 
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1.3 
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201(134) 
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1.5 
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Table 6: Psedoscalar fp decay constants for D mesons in (MeV) 



fp{D) fp{Ds) fp{D,)/fp{D) 





Our 


Others 


Our 


Others 





ur 




Others 


CPP^=0.1 


154(120) 




169(131) 




1.097( 


1.0^ 


U) 




0.3 


197(155) 


230 [59J 
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0.5 
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0.7 
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235 [60J 


1.098 
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M) 


1.15[60] 


0.8 


257(202) 


208 [6T] 
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241 [61J 


1.098 


1.0^ 


M) 


1.164[6T] 


0.9 


265(208) 


201^62] 


291(226) 


249162J 


1.098 


;i.o^ 


^5) 
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;i.o^ 


^7) 


1.23[65J 


1.5 


297(233) 




327(254) 




1.101( 


1.0^ 


^8) 





and following the fact that any c-number, a, commutes with the Hamiltonian, i.e. aH^ 
we express, 



RnP/v{Q) = RniO) 



1 + iSF)p/v 



{Mnv - Mnp) 



n,CW 



(14) 



Here {SF)p = — | and {SF)y = ^ are the spin factor corresponding to the pseudoscalar (J = 0) 
spin coupling and vector (J = 1) spin coupling respectively [ID]. Mn,cw and -Rn(O) are spin 
average mass and the normalized spin independent wave function at the origin of the meson 
state respectively. It can easily be seen that this expression given by Eqn [H] is consistent with 
the relation 

3Rv{0) + Rp{0) 



m 



(15) 



given by [HI \5U\ for nS states. The decay constants by incorporating first order QCD correction 
to the Van Royen-Weiskopff formula are given by [571 EH] ? 



fp/v{nS) 



R 



■nP/V 



(0) 



■nMnp/v 

where, the first order QCD correction factor, C{as) is expressed for the Qq system as 



(16) 



C{as 



1 + 



TT 



mi 



m2 , mi 
m- 



:V,P 



mi + m2 m2 



rria 



17) 
We 



Here 5^ = | [53, [66] and 5^ = 2 [571 ISHl [66] • For the Qq system, mi = mq and m2 
re-examine the predictions of the decay constants fp and /y under different potential schemes 
(by the choices of different v) with and without the QCD correction. Our computed results up 
to 35* states of the D and Dg systems are tabulated in Tables O The ratio of ^fp^£) for 1>S' state 
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Table 7: The root mean square radii (in fm) of D and Dg systems. 
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2P 
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1, 
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2, 
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1, 
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3, 
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2, 
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CPP, 
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.1 


1, 


.255 


3 


.749 
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2 


.819 


5.780 


4 


.381 


5. 


.914 




0, 


.3 


0, 


.979 


2, 


.532 


4.182 


1 


.855 


3.500 


2 


.646 


3, 


.395 




0, 


.5 


0, 


.848 


2, 


.015 


3.161 


1 


.479 


2.646 


2 


.023 


2, 


.524 




0, 


.7 


0, 


.769 


1, 


.710 


2.581 


1 


.267 


2.172 


1 


.682 


2, 


.058 




0, 


.8 


0, 


.740 


1, 


.599 


2.374 


1 


.192 


2.005 


1 


.564 


1, 


.897 




0, 


.9 


0, 


.715 


L 


.506 


2.203 


1 


.130 


1.867 


1 


.466 


1, 


.766 




1, 


.0 


0, 


.694 


1, 


.427 


2.058 


1, 


.077 


1.751 


1 


.385 


1, 


.656 




1, 


.1 


0, 


.676 


1, 


.359 


1.935 


1 


.033 


1.652 


1 


.315 


1, 


.564 




1, 


.3 


0, 


.646 


1, 


.247 


1.735 





.960 


1.493 


1 


.205 


1, 


.417 




1, 


.5 


0, 


.623 


1, 


.159 


1.581 





.905 


1.371 


1 


.120 


1, 


.305 



is tabulated against different choices of u in Table [61 The present results are in accordance with 
other predictions as seen from the the pseudoscalar decay constant fc and fn^. 

4 Root mean square Radii of the D and Dg meson states 
and Average quark Velocities 

The mean square size of the meson is an important parameter in the estimations of hadronic 
transition widths [671 El [Z2], while the average velocity of the quarks within a quark- ant iquark 
bound state is important for the estimation of the relativistic corrections and useful particularly 
in the NRQCD formalism. It is also important in the estimation of their production rates [73j . 
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Table 8: Mean square velocity of the quark within D and Ds states. 



State Model IS 2S 3S IP 2P ID IF 



0, 


.1 


0. 


,197 


0. 


,091 


0, 


.066 


0, 


.098 


0, 


.069 


0, 


.077 


0. 


,069 


0, 


.3 


0. 


,307 


0. 


,208 


0, 


.184 


0, 


.220 


0, 


.190 


0, 


.208 


0. 


,208 


0, 


.5 


0. 


,397 


0. 


,334 


0, 


.331 


0, 


.341 


0, 


.332 


0, 


.353 


0. 


,372 


0, 


.7 


0. 


,472 


0. 


,469 


0, 


.504 


0, 


.459 


0, 


.492 


0, 


.506 


0. 


,556 


0, 


.8 


0. 


,506 


0. 


,539 


0, 


.599 


0, 


.517 


0, 


.578 


0, 


.584 


0. 


,653 


0, 


.9 


0. 


,537 


0. 


,610 


0, 
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0, 
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0, 


.666 


0, 
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0. 
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1 


.0 


0. 


,566 


0. 
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0, 
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0, 
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0, 
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0, 


.740 


0. 


,853 


1 


.1 


0. 


,593 


0. 


,754 


0, 


.914 


0, 
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0, 
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0, 


.818 


0. 
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1, 


.3 


0. 


,642 


0. 
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1, 
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0, 
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1, 


,041 


0, 
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1. 
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1, 


.5 


0. 


,685 


1. 
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1, 


.381 


0, 
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1, 
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1, 
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1. 
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.1 


0. 
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0. 


,060 
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.065 


0, 
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0. 
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0, 


.3 


0. 


,205 


0. 


,136 


0, 
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0, 


.145 


0, 


.124 


0, 
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0. 


,135 


0, 


.5 


0. 


,263 


0. 


,219 


0, 


.217 


0, 


.224 


0, 


.217 


0, 


.231 


0. 


,243 


0, 


.7 


0. 


,313 


0. 
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0, 


.330 


0, 


.302 


0, 


.323 
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0. 


,364 
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0. 


,335 


0. 
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0, 
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0, 
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,428 


0, 


.9 


0. 
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0. 
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0. 


,448 


0, 


.528 


0, 


.413 


0, 


.497 


0, 


.486 


0. 
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.1 


0. 
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0. 
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0, 


.600 
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.448 


0, 


.559 


0, 


.537 


0. 


,627 


1 


.3 


0. 


,425 


0. 


,592 


0, 


.751 


0, 


.516 


0, 


.684 


0, 


.639 


0. 


,762 


1, 


.5 


0. 


,454 


0. 


,687 


0, 


.909 


0, 


.579 


0, 


.812 


0, 


.737 


0. 


,897 
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Thus with our numerical radial wave functions obtained for different choices of the potential 
index u, we compute the mean square radii of the meson state as 

{r')^^ = J\'\Rni{r)\'dr (18) 

and the average mean square velocity of the quark/ ant iquark inside the state as [71] 

= T^iEni - (V(r))„0 (19) 

nl ^^^^1 

Here, E^i is the binding energy of the ni^'^ state and {V{r))ni is the expectation value of the 
quark- ant iquark interaction (without spin dependent terms) potential energy in that state. The 
computed root mean square radii and the relative mean square velocities of the bound states 
within the mesons are tabulated in Table [7] and Table [8] with different choices of u respectively. 



5 Inclusive Semileptonic Decay of Open Charm Flavour 
Mesons 

Inclusive widths of the heavy flavor hadrons are examples of the genuine short- distance processes. 
The open charm mesons, decay through, c qi~^i>, where q = d, s. The light d or s daughter 
quark is bound to the initial light quark of the charm meson by the strong interaction to form a 
new hadron X, according to the Feynman diagram of Figl2J 

In semileptonic decays, the two leptons do not feel the strong interaction, and are thus free of 
strong binding effects. Therefore, they can be factored out of the hadronic matrix element in the 
amplitude of the semileptonic decay process as 

A = ^K;^7m(1 - 75)/(X|g7"(l - 75)c|Z^), (20) 

where all strong interactions are included in the hadronic matrix element (X|g7^(l — 75)0!/^). 
The amplitude of the semileptonic decay process depends both on the hadronic matrix ele- 
ment and the quark-mixing parameter V^^-the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ment. Thus, the semileptonic charm meson decay process is a good laboratory for both studying 
the quark-mixing mechanism and testing theoretical techniques developed for calculating the 
hadronic matrix element. The hadronic matrix element can be decomposed into several form 
factors according to its Lorentz structure. The form factors are generally controlled by non- 
perturbative dynamics, since perturbative QCD could not be applied directly. 
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Figure 2: Feynman diagram for semileptonic D Decay 



For the present study, we compute the decay width of the D and Dg mesons using the nonper- 
turbative bound state effects. The decay width and branching ratio for the r{D + e"*" + Ve) 

and T{Ds — > + ^+ + z/^) mesons are calculated using the expression given by [Ml [69l 170] . 



1927r3 



cd I 



a. 



f{x) g{x) 



TC 



(21) 



1927r3l'^'=^l 



a. 



f{x) -g{x 



IT 



(22) 



where f{x) 
given by [501 168] 

g(x) 



l — 8x + 8x^ — x^ 



12x^ logo;, and the analytic expression of the function g{x) is 



-15.28a;^ + 48.68a;^ - 60.06a;^ + 35.3x^ - S.llx^ - 1.97a; + 2.41 



(23) 



ml/ {nfJ^Y . Generally, for the calculation of the 



Here, the parameter x is computed as x 
semileptonic decay of the heavy flavour mesons, the is taken as the model mass parameters 
coming from the fitting of its mass spectrum. However, taking into account of the binding energy 
effects of the decaying heavy quark within the potential confinement scheme, we consider the 
decaying heavy quark mass as the effective mass of the quarks, rrfj^ . Accordingly, we define the 
effective masses of the quarks in the Qq system as 



m, 



(Ebind) 

rriQ + ruq^ 

{Ebind) 1 

niQ + niq) 



(24) 



to account for its bound state effects. The binding effect has been calculated as {Ebind) = 
MQq — {mQ + mq), where mq and rUq are the model mass parameters employed in its spectroscopic 
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Table 9: The inclusive semileptonic BR of -D — > K'^ + e"*" + and ^ (p + + ui states in %. 



CPP, ^ 


0.1 


0.3 


0.5 


0.7 


0.8 


0.9 


1.0 


1.1 


1.3 


1.5 


Expt.[47J 


BRd 


7.4 


6.7 


6.1 


5.6 


5.4 


5.3 


5.1 


5.0 


4.7 


4.5 


8.6±0.5 


BRd, 


2.77 


2.57 


2.42 


2.29 


2.24 


2.19 


2.15 


2.11 


2.04 


1.98 


2.36±0.26 



study and Mqg is the mass of the mesonic state. The effective mass of the quarks would be 
different from the adhoc choices of the model mass parameters. For example, within the meson 
the mass of the quarks may get modified due to its binding interactions with other quark. Thus, 
the effective mass of the charm quark will be different when it is in cs combinations or in cd 
combinations due to the residual strong interaction effects of the bound systems. 
From the computed inclusive semileptonic decay widths, the Branching ratio of Dq mesons are 
taken here from the relation 

BR = Tsi X r, (25) 

The Lifetime of these mesons {td = 1.04 ps^^ and t^^ = 0.5 ps~^) are obtained from the world 
average value reported by Particle Data Group (PDG-2008)[17j. The computed results of D and 
Ds mesons are listed in Table: [9l Our results are found to be in agreement with experimental 
results at lower potential indexes z/ «i 0.1 to 0.5, which in consistent with the agreement observed 
for their spectroscopy. 

6 Leptonic Decay of the Open Heavy Flavour Mesons 

Charged mesons formed from a quark and anti-quark can decay to a charged lepton pair when 
these objects annihilate via a virtual boson (See FigO]). quark- ant iquark annihilations 




Figure 3: Feynman diagram in standard model for ^ £ u decay. 
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via a virtual to the i'^uli^u) final states occur for the 7T^,K^,Df and mesons. 

There are several reasons for studying the purely leptonic decays of charged mesons [65]. Such 
processes are rare but they have clear experimental signatures due to the presence of a highly 
energetic lepton in the final state. The theoretical predications are very clean due to the absence 
of hadrons in the final state [75]. The total leptonic width of D, Dg mesons are given by 

r(D+ tv,) = ^fljV.,\'mj (l - ^)'m^„ q = d,s (26) 

These transitions are helicity suppressed ; i.e., the amplitude is proportional to m^, the mass of 
the lepton i, in complete analogy to vr^ i'^u. 

The leptonic widths of the charged D and Dg ( I^Sq state) mesons are obtained using Eqnl26] 
employing the predicted values of the pseudoscalar decay constants and /d^ along with the 
masses of the Mu and Md^ obtained from the CPP^ model. The leptonic widths for separate 
lepton channel by the choice of mi=r,fi,e are computed. The branching ratios of the total leptonic 
widths are then calculated using Eqn. [251 The present results as tabulated in Table [10] are in 
accordance with the available experimental values. 



7 Result and Discussions 



The spectroscopic results obtained for open charm {D, Ds) mesons with different choice of the 
confining potential index u from 0.1 to 2.0 are tabulated along with other relativistic quark model 
predictions and with the known experimental states. Our predicted masses of P— wave meson 
state 13P2(2342 - 2514 MeV), l3Pi(2361 - 2542 MeV), l3Po(2335 - 2510 MeV) and liPi(2269 - 
2385 MeV) for the choices of u, 1.0 < z/ < 2.0 are in a accordance with other theoretical model 
prediction [2], [26l [19]. Similar agrement for the predicted masses for the 25", ID, 2P, 3S and 
IP states are also observed in the same range of u values. While the experimental candidate for 
= 2+ P)*(2460), = 1+ P)(2420) [17] and = 0+ state observed in the range 2300 - 2400 
MeV by Belle and Focus [76] lie within our predicted range. In the case of open strange-charm 
mesons {Ds), our predictions within the range of u, 1.5 < < 2.0 for the IP, 25*, ID, 2P, 
35* and IP states of the Dg mesons are in accordance with other theoretical model prediction 
[2l|26lll9]. In particular the experimental states of Ps2(2573) [U] lie within our predicted range 
of l^Pa (2416-2573), P,i (2460) [37] lies in the predicted range of l^P^ (2397-2549 MeV) for 
1.0 < u < 2.0 and the D*^^ (2317) [37] lies close to the predicted range of I^Pq (2317-2350 MeV) 
for 0.8 < z/ < 1.0 of the CPP^ model. The radial excitation of P;(2715) by the Belle grouplTTj is 
found to be close to the predicted 3^5"! state for u = 1.0. Even higher excited states of cs system 
has been observed by the BaBar collaboration [78j with spin parity O"*", l"*" and 2+ etc., with 
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mass at 2856±1.5±5.0 which in our case corresponds to the 2P state with the predicted mass 
range of 23p2(2668 - 2954 MeV), 23Pi(2651 - 2928 MeV), 23Po(2612 - 2858 MeV) and 2iPi(2656 
- 2935 MeV) with the choices of u in the range 1.0 < jy < 1.5. Thus the present study on the 
open charm and open strange-charm mesons using CPP^ model, we have been able to identify 
the recently discovered D— meson states as well as the Dg— meson states. Other predicted high 
angular momentum states ^ > 2 of these mesons are expected to be seen in the future experi- 
ments at BES-III, BaBar, Belle and CLEO collaborations. Our IP— state mass predictions are 
in accordance with the theoretical predictions based on a relative quark model [2] but at higher 
choice of z/ (z/ > 1.5). 

Our results for fp and fv in the potential index ranging from 0.5 to 1.5 are fairly close to the 
known theoretical prediction as seen from Tables [51 The present tabulated results with QCD cor- 
rections (shown in brackets) are in agreement with the experimental values but higher potential 
index beyond u = 1.0. CLEO has reported the first significant measurement of fr)+ = 222. 6± 16.7 
MeV [79] which is close to our predicted value of 227 MeV (without QCD corrections) for u = 0.5 
and that 226 MeV (with QCD corrections) obtained at z/ = 1.3. The accuracy of the previous 
world average has been improved by BaBar with fos = 283 ± 17 MeV [HQ] which within the range 
of values 273 - 283 MeV predicted here for the potential index 0.7 < z/ < 0.8 without QCD correc- 
tion and goes beyond z/ = 1.5 with QCD correction . However the ratio fp{Ds{lS)) / fp{D{lS)) 
is very close to each other between 1.09 to 1.10 (without the QCD correction) and between 
1.082-1.088 (with QCD correction) with changing u from 0.1 to 1.5. The ratio predicted by the 
CPPjy model is thus very close to the ratio predicted by ^59j and ^63j but is lower than the ratio 
of 1.27 as per the recent experimental values of CLEO [79] and BaBar [80] . 
The semileptonic branching ratios of D and Dg mesons computed here using CPPj, model (See 
Table [9]) are all found to be in good agreement with their respective experimental results. It can 
also be seen that the results do not vary appreciably with change in potential index z/, indicating 
lesser influence of strong interaction effects in these decays. Though our predictions for Ds are 
well within the experimental error bar, the branching ratio of D— meson is slightly under esti- 
mated. 

Present study on the leptonic decay branching ratios of D and Ds system presented in Table 
[To] are as per the available experimental limits. The branching ratio in r— lepton channel for 
D and Dg mesons lie within the predicted range for the potential index u ^ 0.3 to 0.5. In the 
case of /i— lepton channel, the experimental value of (4.4 ± 0.7) x 10^^ for D— meson lie in the 
predicted range for the potential index u = 0.3 to 0.5 and that for Dg meson in the potential 
index u = 0.7 to 0.8. Large experimental uncertainty in the electron channel make it difficult for 
any reasonable conclusion. Probably, future high luminosity better statistics and high confidence 
level data sets will be able to provide more light on the spectroscopy and decay properties of 
these open charm mesons. 
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Table 10: The leptonic BR of D and Dg mesons. 





BR^ X 10-3 


BR. X 10-^ 


BRe X 10-8 


D CPP^ =0.1 


1.5 


2.2 


0.5 


0.3 


1.7 


3.6 


0.8 


0.5 


1.6 


4.7 


1.1 


0.7 


1.3 


5.6 


1.3 


0.8 


1.2 


5.9 


1.4 


0.9 


1.0 


6.3 


1.5 


1.0 


0.9 


6.6 


1.5 


1.1 


0.7 


6.9 


1.6 


1.3 


0.5 


7.3 


1.7 


1.5 


0.3 


7.7 


1.8 


Expt.[47j 


< 2.1 


4.4±0.7 






BR^ X 10-2 


BR^ X 10-3 


BRe X 10-^ 


D, CPP^ =0.1 


4.3 


2.5 


0.6 


0.3 


6.3 


4.1 


1.0 


0.5 


7.4 


5.4 


1.3 


0.7 


8.0 


6.4 


1.5 


0.8 


8.2 


6.9 


1.6 


0.9 


8.3 


7.3 


1.7 


1.0 


8.4 


7.7 


1.8 


1.1 


8.4 


8.0 


1.9 


1.3 


8.4 


8.6 


2.0 


1.5 


8.3 


9.1 


2.1 


Expt.lU] 


6.6±0.6 


6.2±0.6 
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